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We present a search for the production of neutral Higgs bosons decaying into t + t~ pairs in pp 
collisions at a center-of-mass energy of 1.96 TeV. The data, corresponding to an integrated luminosity 
of 5.4 fb _1 , were collected by the DO experiment at the Fermilab Tevatron Collider. We set upper 
limits at the 95% C.L. on the product of production cross section and branching ratio for a scalar 
resonance decaying into t + t~ pairs, and we interpret these limits as limits on the production of 
Higgs bosons in the minimal supersymmetric standard model (MSSM) and as constraints in the 
MSSM parameter space. 

PACS numbers: 14.80.Ec,14.80.Da,13.85.Rm 



Supersymmetry (SUSY) [J is one of the extensions of 
the standard model (SM) proposed to address its short- 
comings, such as the hierarchy problem caused by the 
quadratically divergent radiative corrections to the Higgs 
boson mass. In the minimal supersymmetric standard 
model (MSSM), two complex Higgs boson doublets lead 
to five physical Higgs bosons: two neutral CP-even (h, 
H), one neutral CP-odd (A), and two charged Higgs 
bosons (H^). The three neutral Higgs bosons (h,H,A) 



*with visitors from "Augustana College, Sioux Falls, SD, USA, 
''The University of Liverpool, Liverpool, UK, C SLAC, Menlo Park, 
CA, USA, ^University College London, London, UK, e Centro 
de Investigacion en Computation - IPN, Mexico City, Mexico, 
ECFM, Universidad Autonoma de Sinaloa, Culiacan, Mexico, and 
9 Universitat Bern, Bern, Switzerland. 'Deceased. 



are collectively denoted as <f>. At tree level, the Higgs sec- 
tor of the MSSM is fully described by two parameters, 
which are commonly chosen to be the mass of the CP-odd 
Higgs boson, Ma, and the ratio of the vacuum expecta- 
tion values of the two Higgs doublets, tan j3. Radiative 
corrections introduce dependencies on additional MSSM 
parameters. The neutral MSSM Higgs bosons decay into 
t + t~ and 66 pairs with branching ratios of ~ 10% and 
~ 90%, respectively Their production cross section is 
enhanced by a factor that depends on tan /3 with respect 
to the cross section for the SM Higgs boson at the same 
Higgs boson mass. Moreover, for large tan/3, the Higgs 
bosons A and either h or H arc nearly degenerate in mass 
which leads to an effective doubling of a^M^). 

Searches for the production of neutral MSSM Higgs 
bosons have been performed at the CERN e + e~ Col- 
lider (LEP) 0. The CDF and DO Collaborations at 
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the Fermilab Tevatron Collider and the CMS Collab- 
oration at the CERN Large Hadron Collider have ex- 
cluded Ma of up to 300 GeV in a restricted region of 
tan /? ps 30 — 100, by searching for the exclusive processes 
(6)60 -> (6)666 [| and 60 -> 6r+r" 0, and for the in- 
clusive process — > t + t~ 

This Letter presents a search for the inclusive process 
gg, 66 — > — > r + r _ , where the tau lepton pairs are recon- 
structed through their decay into e/i or /xt/i final states, 
and Th represents the hadronic decay modes of the tau 
lepton. The search for r + r~ final states is interpreted 
in a model-independent way and in the context of the 
MSSM. The data were recorded with the DO detector [§] 
at a pp center-of-mass energy of \fs = 1.96 TeV and cor- 
respond to an integrated luminosity of 5.4 fb _1 . This 
represents a significant increase compared to the results 
previously published by the CDF and DO Collaborations, 
which are based on integrated luminosities of 1.8 fb _1 Q 
and 1.0 fb -1 ||, respectively. 

Signal samples arc generated using the pythia [l(| 
Monte Carlo (MC) event generator with the CTEQ6L1 
parton distribution functions (PDF) (llj . Dominant 
background processes comprise Z+jets, TV+jets, and 
multijet production. Background from multijet events 
arises when jets are misidentified as leptons. Additional 
backgrounds include ti and SM diboson production. The 
backgrounds from Z+jets, W+jets, and ti production 
are modeled using ALPGEN [T^], with parton showering 
and hadronization provided by pythia. The ALPGEN- 
generated samples make use of the MLM [3] jet-parton 
matching scheme to improve the jet multiplicity model- 
ing. Diboson processes ( WW, WZ, ZZ) are simulated us- 
ing PYTHIA. In all cases tauola [lj| is used to model the 
tau lepton decays. Simulated events are then processed 
by a GEANT-based 15| simulation of the DO detector, and 
data events from random beam crossings are overlaid to 
model detector noise and multiple pp interactions. Higher 
order quantum chromo dynamics (QCD) calculations of 
cross sections are used to normalize the simulated back- 
ground samples, except for the background from multijet 
production, for which the normalization and differential 
distributions are derived from data. 

Events are selected by requiring at least one single 
muon trigger for the fiTh channel, while for the e\x chan- 
nel, they need to fulfill cither inclusive electron or muon 
trigger conditions. Electrons are reconstructed using 
their characteristic energy deposits, including the trans- 
verse and longitudinal shower profiles in the electromag- 
netic (EM) calorimeter. Muons are identified by combin- 
ing tracks in the central tracking detector with patterns 
of hits in the muon spectrometer. Electrons and muons 
are required to be isolated in the calorimeter and in the 
tracking detectors. 

Tau lepton decays into hadrons are characterized as 
narrow, isolated jets with lower track multiplicity than 
quark or gluon jets. Three types of tau lepton decays 



are distinguished by their detector signature. One-prong 
tau decays consisting of energy deposited primarily in 
the hadronic calorimeter associated with a single track 
(tc i'-like) are denoted as r-type 1; r-type 2 corresponds 
to one-prong tau decays with energy deposited in both 
the hadronic and EM calorimeters, associated with a sin- 
gle track (p^ v- like); and r-type 3 are multi-prong decays 
with energy in the calorimeter and two or more associated 
tracks with invariant mass below 1.7 GeV. A calibration 
for the energy of Th candidates measured in the calorime- 
ter is derived from data. It is based on the ratio of the 
calorimeter energy and the transverse momentum, pr, 
measured in the tracking detector for the Th candidates. 
The ratio is adjusted in the simulation to match the data 
as a function of the fraction of the Th energy deposited 
in the EM calorimeter. 

A set of neural networks, one for each r-type, is ap- 
plied to discriminate hadronic tau decays from jets [16( ■ 
The input variables are related to isolation and shower 
shapes, and exploit correlations between calorimeter en- 
ergy deposits and tracks. When requiring the neural net- 
work discriminants (NN T ) to be JVJV r > 0.9 for r-types 1, 
2 and JVJV r > 0.95 for r-type 3, approximately 67% of 
Z/-f* —> t + t~ events are retained in data, while 98% of 
the multijet background events are rejected. 

A series of selections is used to reduce the background 
from Z+jets, V^+jets, and multijet production. The 
Z/7* — > t~t + process differs from a Higgs boson signal 
only through the mass and spin of the produced reso- 
nance and cannot be further reduced. One isolated muon 
with p^ > 15 GeV and an isolated hadronic tau lepton 
with transverse energy Ej, > 12.5 GeV (r-types 1,2) or 
EJp > 15 GeV (r-type 3) are required in the \iTh chan- 
nel. The muon and the Th must be oppositely charged, 
where the charge of the Th candidate is determined by 
the curvature of the associated track. For r-type 3 the 
charge is obtained by summing over all tracks associated 
with the Th,. The pseudorapidity 77 [l7| is required to be 
\r\fj\ < 1.6 for muons and |r/ T | < 2.5 for tau leptons. The 
transverse momentum sums of all tracks associated with 
the Th candidate, pj,, are required to be greater than 
7, 5, 10 GeV for r-types 1, 2, and 3, respectively. At 
least one hit in the active layers of the DO silicon ver- 
tex detector is required for the tracks associated with 
the Th- The Th and the muon are required to originate 
from the same pp vertex and must be separated from 
each other by ATI = ^(Ar/) 2 + (A^) 2 > 0.5, where 
Aip is the difference in azimuthal angle. This require- 
ment suppresses the Z/j* — > background. The 
transverse W boson mass in W — > Iv events is given 



by M|? = yJ2p l T E T [l~coa(A(p(£,E T )} with I = e,/x. 
The components fS x and $ y of the missing transverse en- 
ergy, are computed from calorimeter cells and the 
momenta of muons, and corrected for the energy re- 
sponse of electrons, tau leptons, and jets. We require 



5 



M^ 1 < 50 GeV to reject VF(— s- /xi^+jets events where 
jets are misidcntificd as Th candidates. 

In the e/x channel, events with at least one muon with 
p T > 10 GcV and |^| < 1.6, and an oppositely charged 
electron with p T > 12 GeV and \-q e \ < 2 are selected. 
The e/j, pair formed by the leptons with the highest pr 
are selected as a candidate; they must be separated by 
A1Z > 0.4. To reject Z — >• ^7 events, an electron 
candidate is rejected if it shares the same track with a 
muon. Multijct background and W boson production 
are suppressed by requiring the mass of the e\i pair to be 
larger than 20 GcV and $ T + P T + Pt > 65 GcV - Back- 
ground from VF+jets production is reduced by requiring 
mm{M T v ,M^} < 10 GcV. The difference in the az- 
imuthal angle, A<p(£,$t), has to be < 0.3 where I = e, \i 
is the lepton with the smaller px- This requirement re- 
jects background from WW, ti, and W+jets production. 
Requiring the scalar sum of the transverse momenta of all 
jets to be < 70 GeV rejects a large fraction of ti events. 

To determine the expected background contribution 
from multijet production in the fiTh channel, two NN T 
regions are selected in addition to the high JVJV T "sig- 
nal" region defined previously: the "medium" region in 
the range 0.25 < JVJV T < 0.75 and the "low" region with 
NN T < 0.1. The samples are further divided depending 
on whether the muon and the Th candidate have the same 
or opposite charge. Background from VF+jets production 
in these samples is reduced by requiring M^ u < 50 GeV. 
The transverse mass is calculated from the $t an d from 
the azimuthal angle A<p(/j,,$t) between the direction of 
the muon transverse momentum p T and the $t- The 
estimated contribution from MC-simulated background 
processes is then subtracted from the resulting distribu- 
tions, and the shape of the multijct background is de- 
rived from the distributions of same-sign /it^ pairs with 
JVJV T > 0.9. Multijet events mainly populate the low JVJV T 
rcgion, and the ratio of opposite to same-sign [iTh pair 
events in this region is extrapolated to yield the normal- 
ization of multijet events in the signal sample. This esti- 
mate of the multijet background contribution is verified 
by an independent method which uses the medium JVJV r 
region. The difference between the estimates obtained 
by the two methods is used as systematic uncertainty on 
the multijet background. 

Multijet background in the e/i channel is determined 
by applying the same selection criteria as for signal apart 
from the electron likelihood and muon isolation criteria, 
which arc inverted. The normalization is then taken from 
the ratio of the numbers of events in the opposite and 
same-sign samples. 

We search for an enhancement from a t + t~ resonance 
above the expected background in the distribution of the 
visible mass M vis = y/ (P Tl + P T2 + Pt) 2 , which is calcu- 
lated using the four- vectors of the measured tau lepton 
decay products, P T12 , and the missing transverse mo- 
mentum, Pt = (JPt ■, $x j $y j 0) . In the e/i final state, the 
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FIG. 1: Distributions of M v i s in the (a) /iTh and (b) efi chan- 
nels after all selections. The data, shown with statistical un- 
certainties, are compared to the sum of the predicted back- 
grounds for an integrated luminosity of 5.4 fb _1 . The Higgs 
boson signal for M</> = 120 GeV is normalized to a production 
cross section of cr^, = 50 pb. All entries exceeding the range 
of a histogram are added to the last bin. 



four- vectors P T1 2 are calculated using the reconstructed 
electron and muon, respectively. After imposing all se- 
lection requirements, the M v ; s distributions for the /xt>, 
and eyu final states are shown in Fig. [T] Table fl] gives the 
yields of the predicted background and of data, summed 
over the M v { s distributions shown in Fig. [TJ 

Several sources of systematic uncertainty affect both 
the signal efficiency and background estimation. Both 
uncertainties that modify only the normalization and un- 
certainties that change the shape of the A/ V i S distribution 
are taken into account. Those that affect the normaliza- 
tion include the integrated luminosity (6.1%), muon iden- 
tification efficiency (2.9%), r/j identification (12%, 4.2%, 
and 7% for r-types 1, 2, and 3, respectively), efficiency 
to reconstruct the track (1.4%), electron identifica- 
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TABLE I: Expected number of events for backgrounds, num- 
ber of events observed in data and efficiency, relative to all r 
lepton decays, for a signal with M$ = f 20 GeV summed over 
the M v i s distributions shown in Fig.[T] The total uncertainties 
are also given. 



Channel 




e/i 




69f4±59f 


697 ± 55 


Multijet 


972 ± 98 


53 ±8 


W — > ev,fiu, tv 


363 ± 60 


f9±5 


Z/j* — > e + e~ , n + n~ 


353 ± 32 


34 ±6 


Diboson + it 


f 80 ± f 2 


27 ±5 


Total Background 


8782 ± 603 


830 ± 56 


Data 


8574 


825 


Efficiency (%) 


f.f6±0.03 


0.20±0.0f 



tion efficiency (3.5%), PDF uncertainty on the accep- 
tance (4.6%), the uncertainty on the Z+jets cross sec- 
tions (5%), the lU+jets normalization (10% for en and 
20% for fiTh), tt cross section (10%), diboson cross section 
(6%), muon and electron trigger efficiencies (both 5%), 
jet energy scale (1.5%— 2%), and the modeling of the mul- 
tijet background (9.1%, 17.7%, and 12.5% for r-types 1, 
2, and 3, respectively). Uncertainties arising from mod- 
eling of the Z boson transverse momentum and the 
energy scales modify the shape of the M v - ls distribution. 
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Expected 

Observed 
Expected ±1 s.d. 
Expected ±2 s.d. 



O 10" 1 
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WUGeV) 



FIG. 2: Model-independent expected and observed 95% C.L. 
upper limits on the product of production cross section and 
branching ratio for inclusive pp — s> <f> — s> t + t~ production as 
a function of Ms, assuming a SM total width for the Higgs 
boson. The ±1,2 standard deviation (s.d.) variations of the 
expected limits are shown as bands. 



The M v j s distribution is used to calculate upper lim- 
its on the cross section based on a modified frequentist 
method with a log-likelihood ratio test statistics [1S|] and 
a profiling technique to reduce the impact of system- 
atic uncertainties [19|. The value of CL S , is calculated 



TABLE II: Upper limits on the expected and observed cross 
sections (in pb) multiplied by the branching ratio for <j> — > 
t + t~ at the 95% C.L. as a function of Ms (in GeV). 



Ms 


Observed 


-1 s.d. 


Expected 


+ 1 s.d 


90 


14.7 


13.8 


19.2 


27.1 


100 


14.4 


7.00 


10.1 


14.0 


120 


5.22 


2.58 


3.53 


5.01 


140 


2.06 


1.14 


1.60 


2.23 


160 


1.23 


0.75 


1.07 


1.50 


180 


0.80 


0.50 


0.73 


1.01 


200 


0.55 


0.39 


0.54 


0.76 


220 


0.40 


0.33 


0.45 


0.64 


240 


0.36 


0.26 


0.37 


0.53 


260 


0.29 


0.23 


0.32 


0.45 


280 


0.23 


0.19 


0.27 


0.38 


300 


0.23 


0.18 


0.25 


0.36 



as CL S = CL s+ b/CL(,, where CL s+ f, and CI* are the 
p- values under signal+background and background-only 
hypotheses, respectively. The expected and observed lim- 
its arc calculated by scaling the signal until 1 — CL S 
reaches 0.95. The combined limits on the product of pro- 
duction cross section and branching ratio into tau lepton 
pairs are given in Fig. [2] and Table ITll as a function of M^. 
The combined limits assume a scalar resonance with the 
decay width of a SM Higgs boson, which is negligible 
compared to the experimental resolution on M v - ls . 

In addition to Ma and tan /3, the masses and couplings 
of the Higgs bosons in the MSSM depend on additional 
parameters through radiative corrections. The produc- 
tion cross section limits are therefore translated into ex- 
clusions in the tan/3 versus Ma plane for two represen- 
tative MSSM scenarios assuming a CP-conserving Higgs 



sector 12011 . the m™ ax scenario [2l| and the no-mixing sce- 
nario [22| with a Higgs mass parameter fj, = +200 GeV. 
The signal cross sections and branching ratios are cal- 
culated using the feynhiggs [23[ program, where the 
gg — > 4> production cross section is taken from 24[ and 
the bb — > 4> production cross section from [ifl . 

At large values of tan/3, the Higgs boson width in- 
creases with tan /3 and can become significantly larger 
than the value in the SM. This effect was previously 
studied by convolving a relativistic Brcit-Wigncr func- 
tion with the next-to-leading order calculation of the sig- 
nal cross section from FEYNHIGGS as a function of M^ 
and tan /3 Q. In the {Ma, tan/3) region where this anal- 
ysis sets 95% C.L. limits, and for yu = +200 GeV, the 
Higgs boson width is smaller than O.IM^, and less than 
half of the experimental resolution on A/ V i S . The sig- 
nal cross section in this channel is largely insensitive to 
sign(/i). The ratio of the gg —> <fi and bb — > <f> cross sec- 
tions also depends on tan /?. For this inclusive search, the 
difference between the efficiencies of the two production 
mechanisms is small and can be neglected. 

The region in the MSSM parameter space excluded at 
the 95% C.L. is shown in Fig.[3]up to M A = 300 GeV. For 
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100 150 200 250 300 100 150 200 250 300 
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FIG. 3: Expected and observed exclusion regions at the 95% C.L. in the plane of tan /3 versus Ma for the (a) m™ ax and (b 
no-mixing scenarios with fi = +200 GeV. The regions excluded by the LEP Collaborations and the CMS Collaboration @ 
are also shown. 



Ma ~ 140 GeV, the expected exclusion reaches tan (3 ss 
30, which is comparable to recent limits obtained in 
The upper limits on the product of the pp production 
cross section for a neutral Higgs boson and branching 
ratio into tau leptons represent the most stringent limits 
to date. 
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